Optimization of field-dependent nonperturbative renormalization group flows 
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We investigate the influence of the momentum cutoff function on the field-dependent nonpertur- 
bative renormaUzation group flows for the three-dimensional Ising model, up to the second order 
of the derivative expansion. We show that, even when dealing with the full functional dependence 
of the renormalization functions, the accuracy of the critical exponents can be simply optimized, 
through the principle of minimal sensitivity, which yields f — 0.628 and 77 = 0.044. 
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Many phenomena in high-energy as well as in statis- 
tical physics cannot be tackled using perturbative meth- 
ods, either because they come under the strong coupling 
regime of their field theory or because they appear gen- 
uinely nonperturbative like, for instance, confinement in 
QCD or phase transitions induced by topological defects. 
The Wilsonian renormalization group concept, resting 
on progressive integration of high-energy fluctuations , 
has rooted the nonperturbative renormalization group 
(NPRG) formalism. The latter has arousen a renewed 
attention, since the idea emerged 0, 0i to apply this 
procedure to the Legendre transform T of the free energy 
— the Gibbs free energy — rather than to the Hamilto- 
nian. This method has allowed to broach notoriously dif- 
ficult problems, such as the abelian Higgs model relevant 
for superconductivity 0, the phase diagram of Hea 0, 
the Gross- Neveu model in three dimensions 0, H L frus- 
trated magnetsjQj or reaction-diffusion processes 3, ^| 
(see 01 and |l3l| for reviews). 

The NPRG formalism consists in building a sequence 
of running effective average actions Tk , which are infrared 
regulated by a momentum cutoff and thus only include 
fluctuations with momenta larger than the running scale 
k. Tk continuously interpolates between the microscopic 
action Tk=A = S and the Gibbs free energy ^k=cj = T. 
Its flow with k is governed by an exact equation |l4| : 



(2) 

where F). is the full inverse propagator and Rk is the 
infrared cutoff function .38J . 

Though Eq. is exact, hence preserving all non- 
perturbative features of the model, it is functional and, 
obviously, cannot be solved exactly. Any practical calcu- 
lation requires to truncate Tk, most commonly through a 
derivative expansion [l^ . Furthermore, when the num- 
ber of invariants or of renormalization functions grows, 
one often resorts to an additional truncation to simplify 
the numerical task, which consists in field expanding the 
different renormalization functions ( |l4l | and see below). 

Naturally, these expansions raise questions as for their 
convergence and accuracy, since these properties entirely 
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condition the reliability of the method , an d rnany works 
have thus been devoted to their study [11 H El H E 
Eol I21L Esj . These issues have appeared intimately 
related to that of the influence of the cutoff function Rk ■ 
Although the exact solution F = limfc^o Tfc does not de- 
pend on Rk, any truncation breaks this invariance. This 
suggests that the choice of the cutoff may be optimized, 
and various criteria have been proposed, such as the op- 
timization of the rapidity of the convergence of the field 
expansion flil ^3 . ^ t he ^ maximization of the "gap" of the 
propagator 20, 25, 26, or the principle of minimal 
sensitivity (PMS) JS, 29,^0). AU these analyses rely on 
the field expansion of the renormalization functions. 

However, a field expansion cannot always be per- 
formed. First, it is meaningless in dimensions where 
the field canonical dimension vanishes, since all field 
powers then become equally relevant. More importantly, 
some phenomena require, in principle, a functional de- 
scription, for instance if the effective potential develops 
non-analycitics, which can occur for some disordered 
systems |31, 32, 33] and continuous growth models [3^ . 
Since the optimization is crucial to control approxima- 
tions, one has to dispose of an efficient procedure to 
deal with the field-dependent renormalization functions. 
However, this had so far never been investigated as it 
involves a sizeable numerical task. In this article, we 
show that, even without truncating the field-dependence, 
the PMS allows to simply optimize the choice of Rk, and 
we provide the corresponding optimal critical exponents 
for the three-dimensional Ising model, up to the second 
order of the derivative expansion. 

Let us first draw the framework of this calculation. The 
derivative expansion physically rests on the assumption 
that the long-distance physics — corresponding to the 
low-energy modes (g — > 0) — is well captured by the 
lowest order derivative terms. Hence, for the Ising model, 
the ansatz of Tk at order writes UM : 



Ffc = d 



i'x{Uk{p) + lZkip){V^f}, (2) 



where p = "0^/2 is the Z2 invariant. The running poten- 
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tial Uk{p) describes the physics associated with uniform 
field configurations while the field renormalization func- 
tion Zk{p) renders the effect of slowly varying fields. The 
flow equations of the running functions U]^ and Zj. follow 
from Eq. J^l, respectively through its evaluation at a 
uniform field configuration and through isolating the 
dependence of its second functional derivative with re- 
spect to '!/'('?) Etnd 4'{~q) T he running anomalous 
dimension is then defined by |12l |: 

rjk = -dt InZk.o, (3) 

where t =^ ln(A;/A) and Zk.o = Zk{po) is the field renor- 
malization at the running minimum po(fc) of the poten- 
tial. We introduce the dimcnsionless and rcnormalized 
variables p = Zkfl p, Ufc = k^'^ Uk and Zk = Z^j^ Zk, 
which are more convenient for the search of fixed points. 

The field expansion would then consist in expanding 
Uk and Zk in powers of the invariant p around the run- 
ning minimum po^k) of the potential — the latter point 
conveying nice convergence properties (lil |20| — which 
would write for a generic function hk'- 
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We here do not proceed to such an expansion. This 
expression is given for completeness since we shall com- 
pare, in the following, our results with those obtained in 
|29j ensuing from a field truncation of the renormaliza- 
tion functions Uk and Zk to the tenth and ninth order 
respectively — that is p„ — 10 and = 9 in Eq. 

To determine the critical exponents of the Ising model, 
we numerically integrate the flow equations dtUk{p) and 
dtZkip) from an initial microscopic scale A — at which 
the potential is quartic ua{p) = A/2 {p—po)^ and the field 
renormalization is constant za{p) = 1 — to the physical 
scale fc — > 0. For large values of the initial parameter 
po(A), the running minimum po{k) of the potential flows 
to infinity such that the dimensionful renormalized min- 
imum — the magnetization — acquires a finite positive 
value (broken phase). For small initial po{A), the running 
minimum po{k), and thus the magnetization, vanishes at 
a finite scale (symmetric phase). In between, there 
exists a critical initial parameter Pn'^(A) for which the 
running potential reaches a fixed point [35j . The critical 
exponent rj is then given by the fixed point value rj* of 
the running anomalous dimension ijk l35l | . The criti- 
cal exponent v — describing the divergence of the corre- 
lation length near criticality — is related to the negative 
eigenvalue of the stability matrix, corresponding to the 
linearized flows of the renormalization functions on the 
field mesh, in the vicinity of the fixed point. 

Setting Zk{p) = 1 for all k, i.e. neglecting the field 
renormalization constitutes the Local Potential Approx- 
imation (LPA). Within the LPA, the anomalous dimen- 
sion T] remains zero. Including a running field renormal- 
ization coefficient Zk.o, independent of the field, allows 



one to refine the LPA by providing a non-trivial — al- 
though rough — determination of rj. This approximation 
is referred to as First Order Approximation (FOA) in the 
following. The next step then consists in incorporating 
the full field-dependence of the renormalization function 
Zk{p), which is here called Second Order Approximation 
(and denoted SOA). 

We successively study these three levels of approxi- 
mations, and analyze for each the influence of Rk. We 
consider an exponential cutoff — which achieves an effi- 
cient separation of the low- and high-energy modes [l^ 
— parametrized by a free amplitude a: 
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where r{y) = Rk{q^)/{Zkfi q^) and y = q^/k"^ 
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FIG. 1: Variations of v with a at the LPA. The dashed line 
represents the results obtained with a tenth order truncated 
potential j2fll |. dots those obtained keeping the full functional 
dependence of Uk{p) (present work). 



For each parameter a, the initial value of po{A) is fine- 
tuned to reach the critical regime and the associated ex- 
ponents are computed. The optimal a is then determined 
through the PMS, that is for a specific exponent x, the 
parameter apMs for which x is stationary is selected: 



dx{a) 



da 



= 0, 



(6) 



and Xpms = x(q^pms) is identified as the optimized expo- 
nent. The hope is that imposing a constraint satisfied 
by X in the exact theory improves the approximate 
determination of this quantity. The obvious drawback 
of this criterion is that Eq. may exhibit many 
solutions. But in fact, this can be easily circumvented 
by complementing the PMS with a few basic rules that 
allow to simply discriminate between multiple solutions 



We can now set out the results, which are collected in 
Table m along with the best theoretical determinations. 
Figure n displays the variations of with the parameter 
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a within the LPA, together with the analogous curve 
ensuing from the field truncation of the potential (at 
p„ = 10) 29] . Both v(q) curves — with and without 
field truncation — appear to coincide within a few tenths 
of percents. This agreement confirms the convergence of 
the field expansion and thus validates the former analy- 
sis |2^. Moreover, the function v(a) exhibits a unique 
extremum, which hence embodies the stationary solu- 
tion of the PMS. The corresponding optimal exponent 
is i^pMs — 0.6503. 

Let us emphasize that within the LPA, v{a) always 
overvalues the "expected" exponent v = 0.6304(13) (ac- 
cording to the best theoretical estimates) , and that z^pms 
corresponds to the minimum. Thus, as already high- 
lighted in 29], the PMS selects the most accurate ex- 
ponent achievable within a given approximation and is 
therefore equivalent to optimizing the precision. 
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full truncated 


full truncated 


LPA 
(a) FOA 
SOA 


0.6503 0.6506 
0.6291 0.6267 
0.6277 0.6281 



0.1058 0.1058 
0.0443 0.0443 


(b) 7-loop 

(c) MC 


0.6304(13) 
0.6297(5) 


0.0335(25) 
0.0362(8) 



TABLE L Critical exponents of the three-dimensional Ising 
model: a) effective average action method: "full" stands for 
results keeping the full field-dependence of Uk and Zh (present 
work), "truncated" denotes results from field expansions of 
these functions |29|l : V) 6- loop calculations including 7-loop 
corrections |3^ : c) Monte-Carlo simulations [s^- 

We then come to the FOA, for which the anomalous 
dimension becomes nontrivial. The variations of v and 
J] with the parameter a are displayed in Figure [21 com- 
pared with the curves following from a truncated poten- 
tial (at j)u — 10) with a field renormalization coefficient 
■^fc,o (which corresponds to = in Eq. |0J) 2.9] • Both 
approaches yield very similar results. Let us note that, 
for both, the anomalous dimension is defined by the fixed 
point value of ry^ evaluated at the minimum of the po- 
tential (see Eq. (j2Jl). Hence the corresponding results 
are not expected to differ much. Indeed, Figure |21 shows 
that 7?(a) computed from both procedures coincide. On 
the other hand, whether or not a field truncation is im- 
plemented alters the determination of v. Only values of 
the potential and its derivatives at the minimum enter 
the computation of v when the latter is field expanded, 
whereas, when integrating the field-dependent flows, the 
determination of this exponent mixes values on the whole 
field mesh considered and thus incorporates a richer in- 
formation. However, the v(a) curves from both proce- 
dures also reveal a good agreement, below the percent 
level (see Figure |21 and Table QJ. Hence the smallness of 
the discrepancy seems to indicate that the vicinity of the 




FIG. 2: Variations of v and r\ with a at the FOA. The dashed 
line represents results obtained with a tenth order truncated 
potential and a field renormalization coefficient {2^ , dots 
those obtained keeping the full functional dependence of Uh (fi) 
(present work). 



minimum already captures most of the relevant features. 

We can now discuss the optimization at the FOA. Both 
v(a) and 77(0;) still exhibit a single extremum and thus a 
unique PMS solution: i/pMs = 0.6291 and ?7pms = 0.1058. 
Again, these exponents turn out to minimize the distance 
to the best theoretical estimates (see Tabled, since the 
latters lie below all v{a) and 77(a) respectively and the 
PMS solutions are minima. The same conclusion can 
therefore be reiterated as for the equivalence between the 
PMS and the optimization of the accuracy. 

Let us finally investigate the SOA, considering the full 
field-dependence of the field renormalization function 
Zfe(p). This order involves a considerable numerical 
task [s^. Figure shows the variations of the two 
critical exponents with a, together with the analogous 
curves obtained by field expanding both Uk(p) and 
Zkip) (to p„ = 10 and = 9 respectively) [23. First, 
these curves support the previous discussion as for the 
two procedures — with and without truncation — as 
the anomalous dimensions exactly match, and the v 
exponents lie within a percent. Moreover, at this order, 
the i^ia) curvature is reversed such that the unique 
extremum becomes a maximum. However, ^{a) turns 
out to underestimate the "expected" exponent for all a, 
and thus the PMS solution still minimizes the distance 
to the best theoretical value. Indeed, j^pms — 0.6277 
appears in close agreement with the 7-loop estimate 
1/ = 0.6307(13). On the other hand, 77(a) always 
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overvalues the "expected" anomalous dimension and 
Vpms = 0.0443 achieves the minimum. Thus, the PMS 
once again selects the most accurate exponents at a 
given truncation. 
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FIG. 3: Variations of f and r; with a at the SOA. The dashed 
line represents results obtained with a tenth order truncated 
potential and a ninth order truncated field renormalization 
function [29l . dots those obtained keeping the full functional 
dependence of both (p) and (p) (present work) . 



We have proposed an optimization procedure appli- 
cable when dealing with field-dependent NPRG flows. 
The PMS has indeed appeared as an efficient tool, 
allowing to optimize the accuracy of the exponents of 
the three-dimensional Ising model, and thus constitutes 
a useful mean of control of the NPRG formalism even 
when field expansions are not possible or not desirable. 
This procedure would be particularly valuable to con- 
firm the results obtained at the next order (9"*) of the 
derivative expansion [s^l, since verifying explicitely the 
convergence of the field expansion at this order becomes 
extremely tedious. Furthermore, such a procedure would 
ensure one to control the influence of Rk and thus 
to work out reliable results within physical contexts 
requiring a functional description — for instance to deal 
with non-analytical potentials. 
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